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Abstract 

We propose a model for DNA dynamics by introducing the helical structure through twist de- 
formation in analogy with the structure of helimagnet and cholesteric liquid crystal system. The 
dynamics in this case is found to be governed by the completely integrable sine-Gordon equation 
which admits kink-antikink solitons with increased width representing a wide base pair opening 
configuration in DNA. The results show that the helicity introduces a length scale variation and 
thus provides a better representation of the base pair opening in DNA. 
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The B-form DNA double helix molecule is usually modeled by two parallel chains of 
nucleotides known as strands with linkage interms of dipole-dipole interaction along the 
strands and the two strands are coupled to each other through hydrogen bonds between 
the complementary bases ll'l. Molecular excitations in DNA based on the above model 

n n n 

is generally governed by nonlinear evolution equations [2|, y, |J] and in particular by the 
completely integrable sine- Gordon- type equations In the above studies, DNA is 

treated as two coupled linear chains without involving the helical character of its structure. 
However, in nature DNA exists in a double helix form and recently there were attempts 
by few authors to study the dynamics by taking into account the helical character of the 
double helix through different forms of coupling. For instance, Gaeta 3, E Q], Dauxios 
[lo| and Cadoni et al [llj assumed that the torsional coupling between the n*'^ base on 
one strand and the {n + 4)*'* base on the complementary strand is the responsible force 
for the helical nature in DNA and found that the localized excitations are governed by 
.olito. and b.ea.Ke.. . a, QQ and Ca.pa Q howeve. ...oduced ..e Ke„... 

through a proper choice of the coupling between the radial and the angular variables of 
the helix and obtained breathers and kinks. On the other hand, very recently, Takeno [l^ 
introduced helicity in DNA through a helical transformation and obtained non-breathing 
compacton-like modes to represent base pair opening through numerical calculations. 
In this, paper, we propose a model by introducing the helical character in each strand of 
the DNA molecule through a twist deformation of the chain in analogy with the twist in 
cholesteric liquid crystal 16| or orientation of spins in a helimagnet 17| . As an illustration 
in Fig. l(a-c) we have presented a schematic representation of the arrangement of bases, 
spins and molecules respectively in a DNA double helical chain, in a helimagnet and in a 
cholesteric liquid crystal leading to the formation of helical structure. In Fig. 1(a) R and 
R' represent the two complementary strands of the DNA double helix and the dots between 
the arrows represent the hydrogen bonds between the complementary bases. The arrows 
and short lines in Figs. 1(b) and 1(c) respectively represent the spins and molecules at 
different sites and planes in a helimagnet and in a cholesteric liquid crystal. When we go 
along the z-direction, the orientation of spins and molecules are tilted from one plane to 
the next through certain tilt angle. If we join the tips of the arrows representing the spin 
vectors and also, the tips of the molecules they form a helix as shown in Figs. 1(b) and 1(c) 
respectively. 
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FIG. 1: A schematic representation of (a) a DNA double helical chain (b) a helimagnet and (c) a 
cholesteric liquid crystal system 



In a recent paper, one of the present authors studied the nonlinear spin dynamics of 
a helimagnet by incorporating the helicity interms of Frank free energy corresponding to 
the twist deformation which is responsible for helicity in a cholesteric liquid crystal system 
0, The Frank free energy density associated with the twist deformation in a cholesteric 
liquid crystal is given by [p ■ (V x p) — go]^ where the unit vector p represents the director 
axis which corresponds to the average direction of orientation of the liquid crystal molecules, 
go = ^ is the pitch wave vector and q is the pitch of the helix. The discretised form of the 
above twist free energy is written as {[k- (pn x Pn+i)] — ^o}^ where k is the unit vector along 
z-direction. In analogy with the above, we write down the free energy associated with the 
twist deformation in terms of spin vector as {[k- (Sn x Sn+i)] —Qo}^- By taking into account 
the form of free ener gy t he Heisenberg model of Hamiltonian for an anisotropic helimagnetic 



system is written as [17 1 



n 

+/i{[fc-(S„xS„+i)]-go}']. (1) 

In Eq.([T]), Sn = {S!^, S^, S^) represents the spin vector at the n^^ site and the terms propor- 
tional to J and A respectively represent the ferromagnetic spin-spin exchange interaction 
and uniaxial magneto-crystalline anisotropy with the easy axis along z-direction. h de- 
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notes the elastic constant associated with the twist deformation. We identify the above 
hehcal spin chain with one of the strands of the DNA double helical chain. Therefore, 
in a similar fashion we can write down the spin Hamiltonian H2 for another helimagnetic 
system corresponding to the complementary strand with the spin vector Sn replaced by 
S'„. We assume that in the Hamiltonian the exchange, anisotropic and twist coefficients 
as well as the pitch in both the helimagnetic systems are equal. Now, for mapping the 
helimagnetic spin system with the DNA double helical chain we rewrite the Hamiltonian 
by writing the spin vectors as Sn = (S"^, S"^, S'^) = (sin 6^^ cos0„, sin 6'„ sin 0„, cos6'„) and 
Sn = {Sn,S'j^,Sn) = (sin 6*^^ cos0jj, sin^^ sin 0^, cos6'^), where 9n{9'n) and 4>n{4>'n) the 
angles of rotation of spins in the xy and xz-planes respectively. The new Hamiltonian cor- 
responding to Hi is written as 



Hi = [- J{sin On sin On+i cos(0„+i - 0„) 

n 

+ COS On COS On+i } + A COS 6*^ + h{sm On 

X sin On+i sin(0„+i - 0„) - q^Y] , (2) 

We now map the two helical spin systems with the two strands of the DNA double helix 
with the two angles On{0'n) and (j)n{<P'n) representing the angles of rotation of bases in the xz 
and xy-planes of the two strands respectively. A horizontal projection of the n*^ base of DNA 
in the xy and xz-planes is shown in Figs. 2(a) and 2(b). Here Qn and Q'^ denote the tips of 
the n^^ bases attached to the strands R and R' at P„ and respectively. The DNA double 
helix chain is stabilised by stacking of bases through intrastrand dipole-dipole interaction 
and through hydrogen bonds (interstrand interaction) between complementary bases. The 
interstrand base-base interaction or hydrogen bonding energy between the complementary 
bases depends on the distance between them and using the simple geometry in Figs. 2(a,b), 
we can write the distance between the tips of bases as j^, 

{QnQ'nY ~ 2 [sin6'„sin6'^ (cos0„cos0^ + sin0„sin0^) 

— cos On cos 0'n\ ■ (3) 

Now, the above equation represents the hydrogen bonding energy between complementary 
bases and we can write down the Hamiltonian for the interstrand interaction or hydrogen 
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FIG. 2: A horizontal projection of the n''^ base pair in a DNA double helix (a) in the xy-plane 
and (b) in the xz-plane. 



bonds is written as 



Hi2 = f] [sin 9n sin 9'^ (cos < 



I cos 



p'^ + sin (pn sin i 



cos 9n COS 9'^] , 



(4) 



where is a constant. The total Hamiltonian for our helicoidal model of DNA interms of 
the angles of rotation of bases using the above Hamiltonians is written as 



H 



Hi + H2 + H12 



[- J{sin 9n sin 9n+i cos(0„+i - </>„) + cos 9, 
'^sin 9'^^i cos(0 

n+l 4'n) 



n 

X COS 9n^ 



9n+i + sin 9'^ sin 9'^^^ > 
+ cos 9'^ cos 9'^_^i } + h{ [sin9n sin 9n+i 
X sin(0„+i - 0„) - go]^ + [sin 9'^ sin 9'^^^ 
X sin(0^+i - 0^) - go]^} + ^{sin 9n sin 9'„ 
X (cos 0„ COS 0^ + sin 0„ sin 0^) — cos ft^ cos 9'^} 
+A(cos^ 9n + cos^ 9'„ 



Using the equation of motion for the corresponding quasi-spin model 19l| in the limit A » 



(5) 



J, 1], h, we obtain 



2ylcos6'„ and 



2y4'cos6''. Hence, under absolute minima of 
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potential the Hamiltonian ([5]) becomes 



H 



+ 4>'^) + J[2 - cos(0„+i - 0„) 



- cos(0;+i - (f)'J] - r][l - cos(0„ - 
+h{2ql - [sin(0„+i - 0„,) - go]^ 
-[sin(0Ui - 0n) - ?o]^}] , 



(6) 



where / = is the moment of inertia of the bases around the axes at Pn(Pn)- While 
rewriting the Hamiltonian in the above form, we have restricted that the bases are rotating 
in the plane which is normal to the helical axis. In otherwords, we have now restricted our 
problem to a plane base rotator model [6|] by assuming 6 = 6' = ti /2. 



Having formed the Hamiltonian, the dynamics of the DNA double helix molecule can 
be understood by constructing the Hamilton's equations of motion corresponding to the 
Hamiltonian (l6l) as 



l'(j)n = [J + 2/;,cos(0„,+i - </)„)] sin(</)„+i - 
- [J + 2/icos(0„ - (pn-i)] sin(0„ - 
+r]sm{(j)n - (p'n) - 2/igo[cos(0n+i - 
-cos(0„ - 0„_i)], 
= [j + 2/icos(0; 
-[j + 2h 



0n-l) 



n+1 ~ Vn)\ Sini</>„+i 



COS(^0„ - <Pn-l 



sm 0„ - 



+^sin(0^ - 0n) - 2/igo[cos(0;,+i 

-COS(0^ - 0^_i)], 



(7a) 



(7b) 



where overdot represents derivative with respective to time. Eqs. f l7al) and (17bp describe 
the dynamics of the DNA double helix at the discrete level when the helical nature of the 
molecule is represented in the form of a twist-like deformation. 



It is expected that the difference in angular rotation of bases with respect to neighbouring 
bases along the two strands is small j^, Also, very recently Gaeta j^, 9] proposed that 
the helical structure of DNA will introduce qualitative changes only in the small amplitude 
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regime. Hence, under the small angle approximation, in the continuum limit, the discrete 
equations of motion (7a,b) after suitable rescaling of time and redefinition of the parameter 
7] reduce to 



(j)tt = ^^'^^^ (j)zz + r]sm{(f)- (j)'), (8a) 
<P'u = ^-^^^^<PL + VM<P'-<P)- (8b) 
Adding and subtracting Eqs. (IHaj) and (I8bl) and after suitable rescaling of the variable z, we 



obtain 

^H-^.. + sin^ = 0, (9) 

where \l/ = 20 and we have further chosen 2// = —1. Also, while deriving Eq. ([9]), we have 
chosen 0' = —0, because among the possible rotations of bases, rotation of complementary 
bases in opposite direction easily facilitate an open state configuration. Eq. ([9]) is the com- 
pletely integrable sine-Gordon equation which was originally solved for N-soliton solutions 
interms of kink and antikink using the most celebrated Inverse Scattering Transform method 



2l| . For instance, the kink and antikink one soliton solution of the sine-Gordon equation is 
written interms of the original variables as 

(f){z,t) = 2 arctan[ea;p[:' 



where + and — represent the kink and antikink soliton solutions respectively and v is the 
velocity of the soliton. In Fig. 3(a) we have plotted the angular rotation of bases in terms of 
the kink-antikink one soliton solution as given in Eq. ffTOj) by choosing the stacking, helicity. 



moment of inertia and velocity parameters respectively as J = 1.5 eV, h = 3.0 eV, I = 
1.3 X 10~^^(7 cm^ and v = 0.4 cm jiol . The kink-antikink soliton solution which 

can propagate infinite distance and time describes an open state configuration in the DNA 
double helix which is schematically represented in Fig. 3(b). In order to understand the 
effect of helicity on the open state configuration, in Fig. 3(c), we have also plotted the kink 
and antikink one soliton solution of the sine-Gordon equation in the absence of helicity, that 
is by choosing h = (keeping all other parameters values the same), in the solution given 




(b) (d) 

FIG. 3: (a) Kink and antikink one soliton solutions (Eq. (jlOp ) of the sine-Gordon equation when 
helicity is present (/i 7^ 0). (b) A sketch of the formation of open state configuration in terms 
of kink-antikink sohtons in a DNA double helix when helicity is present (h ^ 0). (c) Kink and 
antikink one soliton solutions of the sine-Gordon equation when helicity is absent (Eq. (jlOp when 
h = 0). (d) A sketch of the formation of open state configuration in terms of kink-antikink solitons 
in DNA double helix when helicity is absent (/i = 0). 




FIG. 4: (a) The kink one soliton (Eq. (jlOp ) representing base pair opening at t = 1 for different 
values of helicity. (b) Variation of width of the kink soliton against helicity. 



in Eq. fllOl) . From Figs. 3(a) and 3(c), we observe that when there is helicity in the model 
{h 7^ 0), the kink-antikink soliton is getting broadened. In other words, helicity in DNA 
makes more number of base pairs to participate in the formation of open state configuration 
without introducing any qualitative change in the dynamics. This is also schematically 
represented in Fig. 3(d) which looks evident on comparing Fig. 3(b). In order to highlight 
the above fact, we have separately plotted the kink one soliton solution at a given time 
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{t = 1) for different values of lielicity by choosing h = 0,1, 2, 3, 4 in Fig. 4(a). The increase 
in width against hehcity is exphcitly represented in Fig. 4(b). From the figure it may be 
noted that the increase in hehcity slows down the rotation of bases and makes more and 
more number of base pairs to participate in the open state configuration, thus providing a 
better representation of base pair opening in DNA. Thus, helicity introduces a length scale 
variation in the base pair opening. 

Similar results have also been observed by Dauxios through a perturbation analysis on 
his helicoidal model of DNA and obtained soliton with a much broader width. In order to 
have a more realistic model, dissipative (viscous effect) and noise terms should be added 
to the equations of motion. Experimentally, the life time of soliton in this case is shown 
to be of few nano seconds at room temperature (see for e.g. Ref. 2^). Also, in a recent 



paper, Yakushevich et al 



23| through a numerical analysis, showed that when the viscosity 



is strong the soliton moves a length of only few chain links and it will stop after that. 
On the other hand, when the viscosity is low the soliton passes more than 3000 chain 
links like a heavy Brownian particle which is found to be stable with respect to thermal 
oscillation. When the above two effects are taken into account Eq. ([9]) takes the form 
^« — "^zz + sin^ = e[a;^f + f3r{z,t)] where the terms proportional to a and (3 are related 
to viscous surrounding and thermal forces respectively. A soliton perturbation analysis 2^ 
on the above equation shows that when the viscosity is high the soliton moves for a small 
distance and then stops. But when the viscosity is low the soliton moves for a long distance 
along the chain. The detailed analytical calculations of the above study will be separately 
published elsewhere. 

In summary, we proposed a new helicoidal model to study DNA dynamics by introducing 
the helical character in analogy with the twist deformation in a cholesteric liquid crystal 
system and the spin arrangement in a helimagnet. The nonlinear dynamics of DNA under 
the present helicoidal model is found to be governed by the completely integrable sine- 
Gordon equation in the continuum limit which admits kink and antikink soliton solutions. 
From the nature of solitons, we observe that helicity introduces a length scale variation 
without causing any change in the shape of the soliton. Due to this scaling variation, the 
width of the soliton increases and hence we obtain broader kink-antikinks. In otherwords 
a large number of bases are involved in the base pair opening, thus leading to a better 
representation. This broadened base pair opening may act as a better energetic activator 
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in the case of RN A- Polymerase transport during transcription process in DNA. As the 
continuum hehcoidal model does not introduce qualitative changes in the DNA dynamics, 
we propose to study the full nonlinear dynamics of the helicoidal model of DNA (without 
making small angle approximation) by solving equations fITal) and fl7bl) numerically and the 
results will be published elsewhere. 

The work of M. D and V.V form part of a major DST project. 
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